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l- We apply the projection operator method (POM) to theory and derive both quantum 

• and semiclassical equations of motion for the soft modes. These equations have no time- 

' convolution integral term, in sharp contrast with other well-known results obtained using 

, the influence functional method (IFM) and the closed time path method (OTP). However, 

■ except for the fluctuation force fleld terms, these equations are similar to the corresponding 

equations obtained using IFM with the linear harmonic approximation, which was introduced 
to remove the time-convolution integral. The quantum equation of motion in POM can be 
^ , regarded as a kind of quantum Langevin equation in which the fluctuation force fleld is 

(N . given in terms of the operators of the hard modes. These operators are then replaced with c- 

numbers using a certain procedure to obtain a semiclassical Langevin equation. It is pointed 
(N ■ out that there are significant differences between the fluctuation force fields introduced in 

this paper and those introduced in IFM. The arbitrariness of the definition of the fiuctuation 
'n— «' ' force field in IFM is also discussed. 



§1. Introduction 

In recent years, there have been many studies related to time dependent phenom- 



a^uaJitiim systemsi thetime development of order parameters in phase tran- 
EP'l3)'l§'i)'B)'B'Q)'|)'@)'t3)'llJ the experimental observation of non-exponential 



Oh: 

QJ , ena m 

' sitions,- - -- -- -- - — — ^^^^ ^^^^^^^^^^^^^^^^ ^^^^^ .^.x^x. 

decay in quantum tunneling, the absence of the quantum Zeno efiectin quantum 
^ I field theory, EJ) the time evolution of thfiBosejijiii^ein mndmisate,EJ)'t3) renormal- 

ization in time evolution calculations, ll3)'llZP'll3)' 113) 'E3)'EJ'EHP and so on. In such 
studies, it is often important to evaluate the time evolution within quantum field 
theory. However, the equation of motion in quantum field theory is a nonlinear op- 
erator equation, and hence it is difficult to solve in general. There have been several 
attempts to solve it approximately. 

One attempt to approximately solve the equaiton of motion involves a deriva- 
tion of a semiclassical equation of motion that has a fluctuation force field term, in 
other words, colored noise. This equation can be regarded as a kind of a Langevin 
equation, and it may be possible to incorporate the effects of both quantum and 
thermal fluctuations into it through not only the coefficients in the equation but 
also the fluctuation force field. The Langevin equation is a well-developed tool in 
nonequilibrium statistical mechanics, and there exists many techniques for treat- 
ing it, including mode coupling theory to calculate dynamical critical exponents, c3 
scaling theory of nonequilibrium systems near instability points. Furthermore, 
the effect of the fluctuation force field in the nonlinear Langevin equation is an im- 
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portant topic of study in varioim^fields of physics, e.g. stochastic resonance eEP and 
noise- induced phase transition. E^) Such effects have not been studied in detail in 
cases that quantum field theory is needed to describe the relevant phenomena, for 
example the time evolution of a disoriented chiral condensate, which may be a sig- 
nature of quark-gluon plasmas. The Langevin equation is the starting point in the 
analysis of such phenomena. 

The derivation of Langevin equations in quantum field theory has been carried 
out mainly on the basisx)f ihe iufliience functional method (IFM) and the closed time 
path method (CTP).&|)'H)'i'l§'El)'E|) The resultant Lan gevin equation has two notable 
characteristics. First, it has in general a time-convolution integral term, which is 
often called a memory term. Such a term makes the Langevin equation difficult to 
solve. For thisreason, it is usually removed by employing the quasi-instantaneous 
approximation EJ)'I3^ or the linear harmonic approximation. 13 The meaning of the 
quasi-instantaneous approximation is very clear. BI|) However, it is known that with 
this approximation there can be no dissipation effect (as, for example, in (p^ theory), 
and therefore an additional approximation is employed. Contrastingly, with the 
linear harmonic approximation there is a dissipation effect in cj)'^ theory without the 
use of an additional approximation. However, the validity of this approximation is 
not clear. The second characteristic of the Langevin eqaution obtained using IFM 
and CTP is that the fluctuation force field is introduced using an auxiliary field. It 
is not clear why the auxiliary field introduced in this way can be interpreted as the 
fluctuation force field. £3)' Ej) Furthermore, the definition of the fluctuation force field 
is not unique, as is discussed in this paper. 

The projection operator meth od ( PO M) is ajiother method to derive the Langevin 
equation for quantum systems. o)'I^'I^'E3)'E3)'0'E3''E3) A particular version of this 
method was proposed by Hashitsume, Shibata and ShingiijLi3''EiP and was improved 
by Uchiyama and ShibataE3) and Koide and MaruyamaEB) independently. Using 
the Mori projection as the projection operator, one can derive the Mori equation, 
which is xL,we]I-known linear Langevin equation applicable to classical and quantum 
systems. H'EIP POM is the only method that can be used to systematically derive 
an equation of motion without the time-convolution integral, as far as we know. The 
fluctuation force field at the quantum level is explicitly given as a term represent- 
ing the time variation in the space that is projected out, and there is no ambiguity 
regarding the nature of the fluctuation force field. Therefore, we believe that with 
POM it is possible to derive a new semiclassical Langevin equation that is free from 
the problems mentioned above. 

The purpose of this paper is to derive both quantum and semiclassical Langevin 
equations within quantum field theory using POM and to compare the result so 
obtained with that obtained from IFM. We first derive a quantum Langevin equation 
using POM. In practice, it is difficult to numerically solve such a nonlinear operator 
equation. For this reason, we derive a semiclassical Langevin equation from the 
quantum equation employing a replacement procedure. We apply POM and IFM to 
0^ theory and derive Langevin equations for the soft modes projecting or integrating 
out the hard modes. The method of deriving the semiclassical Langevin equation 
from the quantum Langevin equation is not unique. Therefore, it is useful to compare 
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the semiclassical Langevin equation obtained using IFM with the quantum Langevin 
equation obtained using POM and investigate the differences between their results 
beforehand. In so doing, we show that the result given by POM is similar to that 
given by IFM with the linear harmonic approximation. Furthermore, the differences 
between the fluctuation force fields of POM and IFM can be partially eliminated 
by utilizing the arbitrariness of the definition of the fluctuation force field in IFM. 
Next, we discuss the derivation of the semiclassical Langevin equation in POM. 

This paper is summarized as follows. In §2, POM is reviewed. In §3, we apply 
POM to (j)'^ theory and derive a quantum Langevin equation for the soft modes. 
In §4, the derivation of the semiclassical Langevin equation in IFM is given. This 
section follows the work of Greiner et al.EF The results of §3 and §4 are compared 
in §5. In §6, we derive the semiclassical Langevin equation by replacing the relevant 
operators with c-numbers. The result is compared with the semiclassical Langevin 
equation in IFM. Conclusions and discussion are given in §7. 

§2. Projection operator method 

TheveisioiL of POM that we use in this paper has the following characteris- 
tics.E3)'Ei|)'E|).E6|) First, it can be used in both Schrodinger and Heisenberg pictures; 
in other words, the master equation and the Langevin equation can be treated si- 
multaneously. Second, there is great freedom in the choice of the projection op- 
erators. In fact, this projectiDnonerator method includes both the Mori and the 
Nakajima-Zwanzig methods. Ey^'l^'EHP Third, the equations both with and without 
a time-convolution integral can be treated systematically. Fourth, it can be used 
even when the Hamiltonian is ex p licitly time dependent. In this paper, we derive 
the result in Heisenberg picture. 123 

The starting point is the Heisenberg equation of motion, 

±0{t) = i[H,0{t)] 

= iLO{t), (2-1) 

where L is the Liouville operator. We can rewrite the Heisenberg equation of motion 
by using the projection operator, which has the following general properties: 

p2 = (2-2) 
Q = 1 - P, (2-3) 
PQ = QP = 0. (2-4) 

Then, the Heisenberg equation of motion can be rewritten as 



+Qe-«(-*o)^_J^.LO(to), 



(2-5) 
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where 



U{t,to)= f dTe-'^^*-^'^PiLQe'^'^^^-^\ (2-6) 

Jtn 



I to 

This equation is called the time-convolutionless (TCL) equation because it does not 
include time-convolution integral terms, as is shown below. When the system has a 
dissipation effect, it comes from the second term on the r.h.s. of the equation. The 
third term represents the time variation in the space which is projected out. Note 
that Q operates from the left in this term. We regard this term as the fluctuation 
force term, as in the case of the Mori equation. This equation is still equivalent to 
the Heisenberg equation of motion, and it is difficult to solve without approximation. 
To make an approximation, we have to specify the projection operator. 

We consider the case in which the total system can be divided into two parts, 
the system and the environment. Then, the Hamiltonian is given by 

H = Hs + He + Hi 

= Ho + Hi, (2-7) 

where Hs and He are the unperturbed Hamiltonians of the system and the environ- 
ment, respectively, and Hj is the Hamiltonian that describes the interaction between 
the system and the environment. The self- interactions of the system and the envi- 
ronment are also included in Hj. We assume that the initial density matrix is given 
by the direct product of the system and the environment density matrices: 

P = Ps®PE- (2-8) 

We then define the projection operator P as 

PO = Tie[peO] = {0). (2-9) 

Here, Tr^; is the trace over only the environment degrees of freedom. 

Using the above described projection operator, we can rewrite the TCL equation 
asS 

4o(t) = e*^(*-*o)PiLO(to) 
dt 

i + {L[t, to) — 

^QgiLQ(i-io)__J__i^O(i^). (2-10) 

Here, the new function C(t, to) is defined as 
C{t,to) 



^iLo(t-to) ^-iL{t-to) 



n=l 



to Jto Jto 



1 + E (-^)" / / dt2--- dtnLiih - to)Li{t2 -to)--- Li{tn " to), 



(2-11) 
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where 

LoO = [Ho,0], (2-12) 

LjO = [Hj,0], (2-13) 

Li{t - to) = e'^<)(*-*o)2,^g-iLo(t-to)^ (2-14) 

The operator C{t,to) is a time-ordered function of the Liouville operator. 

For the second term on the r.h.s. of the equation, we expand C{t,to)Q/{l + 
{C{t, to) — 1)Q} up to first order in Lj. Then, we have an approximate representation 
of the TCL equation, 

ft 
Jto 

+QgiLQ(t-to)_J__^^0(^^). (2-15) 

It can be seen that Eq. (2-15) does not contain a time-convolution integral, from the 
form of the full time-evolution operator, e»^(*-*o)^ which operates from the left in the 
second term on the r.h.s. If this equation did contain a time-convolution integral, 
the form of the full time-evolution operator must be e*^(*^*) ^ where s is an integral 
variable. This is the reason that the equation is called the "time-convolutionless 
equation" . 

Now, we expand the third term. It is not clear to what order we should expand 

this term. Here, we apply the fluctuation-dissipation theorem of second kind (2nd f-d 
theorem) to the second and the third terms in order to fix the order of the expansion. 
First, we rewrite the third term as 

= QV{t,to)Q 

1 

^ 1 + {C{t, to) - 1)Q + {V{t, to) - 1) + iC(t, to) - l)Q(V(t, to) - 1) 

(2-16) 

Here, the function V{t, to) is defined as 

V{t,to) = e^Q^Q(f-to) ^-iQLoQ{t-to) 

J^L ft ftl ftn-l ^„ ^„ 

= 1 + ^ / dtl / dt2--- dtnQUf{tn - to)QLY{tn-l - to) 

n=l "^*o "'to Jto 

x---QLf{ti-to), 

(2-17) 
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where 

Lf{t - to) = e^^o(*-*°)L/Qe-^^«(*-*o). (2-18) 

The function io) is an anti-time-ordered function of the Liouville operator. In 
the case of a hnear Langevin equation, hke the Mori equation, the 2nd f-d theorem, 
which relates the coefficient of the equation with the fluctuation force term, is exactly 
satisfied. Therefore, we require it to be satisfied order by order in the perturbative 
expansion. If we carry out the perturbative expansion up to zeroth order in Lj for 
QgiLQ(t-to)_^__^ [that is, if we take C{t,to) = 1 and V{t,to) = 1 in Eq. (2-16)], 
the 2nd f-d theorem is satisfied in the case of the linear Langevin equation, as shown 
in Appendix As a result, we have the approximate TCL equation 

+Qe^^°(*-*«)iLO(to). (2-19) 

Here, we have used the additional condition PLq = LqP for simplicity, because in the 
following calculation, we consider only the case in which this condition is satisfied. 

In this section, we derived the TCL equation from the Heisenberg equation of 
motion. We can derive the TC equation that contains the time-convolution integral 
similarly. (See, for example. Appendix C of Ref. |3^).) In the present paper, we 
consider only the TCL equation because we wish to investigate the validity of the 
approximation employed in IFM in order to eliminate the time-convolution integral. 

We note here that the difference between the expressions "time-convolutionless" 
and "memoryless" must be understood. "Time-convolutionless" does not mean that 
there is no memory effect. In fact, it has been confirmed that the TCL equation has 
a desirable non-Mar kovian effect.©'® Any memory effect that exists in the TC 
equation is also incorporated in the TCL equation through the higher-order terms, 
because these equations are equivalent. 

§3. Quantum Langevin equation in POM 

In this section, we apply POM to (p^ theory. The Lagrangian is given by 

£ = ^d^m - Im^^ - iA<A^ (3-1) 

where the symbol " indicates an operator. We set up the physical picture to be the 
same as that employed in Ref. |5|). In this paper, we treat the case > and do not 
consider the symmetry breaking, to avoid the complication involving the selection 
of the Fock space associated with degenerate vacua. If we were to consider the 
case < 0, the system would have degenerate vacua and would undergo a phase 
transition to the symmetry broken phase. In this case, the expectation value of the 
field operator (p would play the role of an order parameter and would be regarded as 
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a gross variable. For this point, we consider only the > case. We may choose an 
initial conditions for which the order parameter is nonvanishing at t = to, due to its 
history before t = to- The condensate will consist of soft modes. We then introduce 
a cutoff Aj, which is the softest mode among the hard modes. That is, the hard 
modes are defined hy k > Aj. The largest value of k among the modes contained 
in the condensate is assumed to be smaller than Aj. We then impose another set of 
initial conditions in which the hard modes are in thermal equilibrium with a high 
temperature T at t = to- In this situation, we regard the hard modes as irrelevant 
information and carry out a coarse-graining, as in the usual case of Brownian motion. 
Then the hard modes become the origin of dissipation and fluctuation in the soft 
mode equation. 

The field 0(x, t) can be divided into soft modes and hard modes as 



(x,t)= / 'd^k \ K{t)e^'^^ + al{t)e-^'^^] 
Jo V2 27rpcji, 



+ / d3k-=i_[ak(t)e*- + aj,(t)e-*-] 
Jai V2(27r)3wk 

0<(x,t) + 0>(x,t), (3-2) 



where A is an ultraviolet cutoff and u>\^ = + m^. The field (x, t) contains 
only soft modes, and the field (/>> (x, t) hard modes. We consider the case in which 
the hard modes have already thermalized but the soft modes have not yet reached 
thermal equilibrium at t = to- Therefore, the momentum cutoff Aj is the lower 
limit of the momenta that constitute the thermalized hard modes. Then, the hard 
modes play the role of a heat bath for the soft modes. In this case, except for the 
vacuum polarization, which should be renormalized, the temperature T acts as an 
effective cutoff for the hard modes, due to the Bose distribution function. Therefore, 
to ensure that the hard modes have a sufficiently large number of degrees of freedom, 
we must consider the situation that Aj <C TM 

We are interested in only the slow evolution of the soft modes and hence have 
to coarse-grain the fast time dependence of the hard modes. Consider the case of a 
small coupling constant. Then, the time evolution of the soft modes and the hard 
modes is well approximated by that of the free field. The time scale of the slowest 



soft mode is ~ m~^, while that of the slowest hard mode is ~ l/yAj + m^. To 
regard the hard modes as microscopic variables in relation to the soft modes, like 
a kind of random fluctuation, there must be a large difference between the typical 
time scales of the soft modes and the hard modes. For this, the condition m <^ Aj 
must be satisfied. In short, we must impose the following condition for the cutoff 
Ai: 

m^Ai ^T. (3-3) 

In this way, we finally obtain a semiclassical Langevin equation using POM in §6, 
in which the operator-valued fluctuation force terms are replaced with c-number 
stochastic variables. Then, with the aforementioned initial conditions, the Langevin 
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equation describes the time development of the condensate — actuahy, the decay or 
the dissolution of the condensate on a time scale longer than 1/Aj. 
The conjugate field i7(x, is also divided as follows: 



ikxi 



JAi 



2(27r)3 

= n^iy:,t)+n>(yi,t). (3-4) 

The soft modes and the corresponding conjugate fields satisfy the commutation re- 
lations 

[0<(x,t),(^<(x',t)] =0, (3-5) 

[i7<(x,t),i7<(x',t)] =0, (3-6) 

[(^<(x,t),i7<(x',t)] = 4^^)(x-xO, (3-7) 

where 

4^(x-x') = ^/^^'d^ke^'^(— '). (3-8) 

Here, the notation J^^ d^k means that the integration is restricted to the region 

< |k| < Aj. The function S^^' — x') becomes the usual Dirac delta function in 
the Aj ^ oo limit. Hereafter, we call it the "cutoff delta function". 

The Hamiltonians of the system, the environment and the interaction become 



Hs = ^ 


1 


He = ^ 




Hi = ^ 





,2 12 



(3-9) 



r2(x,to) + (v<^>(x,io))' + 

d^x^A[(^^(x, to) + 4|(x, to)^>(x, to) + 6(^|(x, to)<^?>(x, to) 

+^< (x, to)0'^ to) + 4)1 (x, to)] . (3-11) 

Note that the division of H into Hs-, He and Hj is made at the initial time to- 

The initial density matrix is given by the direct product of the system and the 
environment density matrices as in Eq. (2-8). We consider the case in which the 
environment is in thermal equilibrium with the Hamiltonian He^ and thus 

PE = e-P"^/Z, (3-12) 
Z = Tre~^^^, (3-13) 

where /3 = l/T is the inverse temperature. Then, the projection operator is defined 
by Eq. (2-9), with pE given by Eq. (3-12). 
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Substituting the above conditions into Eq. (2-15), we obtain the quantum Langevin 
equation. Here, we summarize the approximations employed in our calculation. 
We prepare the thermal equilibrium state as the initial environment, and expand 
C{t,to)Q/{l + io) — up to first order in Lj. This approximation is made 
for the following reasons. First, it is necessary to take into account at least first 
order in Lj to have a dissipation effect. Second, it was shown in a certain model 
that the system evolves toward the thermal equilibrium state at this level of approx- 
imation. Ej) Furthermore, an approximation similar to that employed in this paper 
is employed in Ref. ^). In that work, the nonlinear Langevin equation is solved nu- 
merically, and the temperature dependences of the order parameter and the effective 
masses of a and vr modes are determined. These temperature dependences are the 
same as those exhibited when the system thermalizes with an initial environment 
temperature T. Third, one of our main purposes is to compare our POM result with 
IFM result, The present approximation is sufficient for this purpose. 

Setting the initial time as to = in Eq. (2-15), we obtain the quantum Langevin 
equation for the operator 0<(x, t), 

(x, t) - A0< (x, t) + (x, t) + ^J (fySf^ (x - y)<^^ (y , t) 



where 



+ 2 J d'y5]^;{^-y){<pUy,t))<P<iy,t) 
+iX^ 1^ ds I d^yid\2Sf^ (x - ys) 

X I ^ ( [4>t> (y 1 , , H> iy2,t)U<{y2,t) 
1 . 

+ ^ ( [(I^L (y 1 , s) (y2 , t)] ) ^< (y2 , s - t, t) 

+^{[^l>iyi,s), 4>ly{y2,m[<fl{yi,s - 1, t), My2,t)]+ 

+ ^ ( [4>o> {yi,s),4>o>{y2,t)])[<fl{yi,s - t,t)Jl{y2,t)]+ 

+^([<^o>(yi,s),0o>(y2,i)]+) 

x['Ao<(yi,s - *),0o<(y2)]['^<(yi,s - t,t),0<(y2,t)]+ 

+I{[4>l>iyi,s),4>l>iy2,t) - {4>l>iy2,m+) 

x['P<{yi,s- t,t),(f><{y2,t)]if<{y2,s - 

= J d^y5f^i^-y){f\{y,t) + f2iy,mo<{y,t) + h4>l<{y,tm-U) 



/i(x,t) = ^0;]>(x,t), (3-15) 
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(a) (b) (c) (d) 





(e) (0 (g) (h) 

Fig. 1 :Diagrams that contribute to the quantum Langevin equation in POM. The 
dashed Hnes represent the soft modes and the sohd hnes the hard modes. 

/2(x,i) = - (0o>(x,O)), (3-16) 

/3(x,t) = ^<^o>(x,t). (3-17) 

Here, the fields with the subscript evolve as free fields, and the newly introduced 
field (^<(x, s,t) is defined as 

Jo ^y2{2^^)'^UJ]^ 

= - — / d^k{(/)<(k, t) cos H iT<(k, i)sina)ks}e''*^'^, 

{2'itY Jo 

(3-18) 

where 



:(k,t) = y'i^(ak(t) + alk(t)), (3-19) 



i7<(k, t) = i^(ak(t) - al^it)). (3-20) 

Furthermore, [ ] and [ ]_|_ in Eq. (3-14) represent the commutation and the anti- 
commutation relations, respectively. The first four terms on the l.h.s. of this equation 
come from the first term on the r.h.s. of Eq. (2T9). The remaining terms on the 
l.h.s. come from the second term, and the terms on the r.h.s. come from the third 
term. The diagrams that correspond to the interaction terms on the r.h.s. of Eq. 
(3-14) are shown in Fig. 1. The diagrams (a) and (b) correspond to the fourth and 
the fifth terms on the l.h.s. of the equation, respectively. The diagrams (c) through 
(h) correspond to the first term through the sixth term within the brackets on the 
l.h.s. of the equation, respectively. 

Now, we call fi the fluctuation force fields. We examine the correlations of 
/j. The expectation values of fi can be calculated using the initial density matrix p. 
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However, it is also possible to replace p with pe, because fi consists of only operators 
of the hard modes. The first-order correlations are 

(/i(x, t)) = (/2(x, t)) = (/3(x, t)) = 0. (3-21) 

The second-order correlations are calculated as follows: 

(/i(xi,ti)/i(x2,t2)) = ^[6Gyi,(xi - X2,ii - ia)^ + 9G^,(xi - X2,ti - i2)Gyi,(0, 0)^], 

(3-22) 

\2 

(/2(xi,il)/2(x2,t2)) = -yGA,(xi-X2,ii-i2)', (3-23) 

iX^ 

(/3(xi,tl)/3(x2,t2)) = -— G,1,(X1 -X2,tl -t2). (3-24) 

The propagator on the r.h.s. of Eqs. (3-22), (3-23) and (3-24) is defined as 

rA I 

GAi{x,t) = i / d^p 3 — {{l + 2n^.p)cosujpt-ismLUpt}e'- 

(3-25) 

It has the property 

GXi^,t) = -GA,{-^,-t). (3-26) 

This propagator includes only the hard modes. The Bose distribution function n^p 
is given by 

The correlation functions in Eqs. (3-22), (3-23) and (3-24) are complex, and 
change form under the exchange xi ^ X2,ti ^ t2- Now, we examine the real 
and the imaginary parts of the correlations. The real parts of the second-order 
correlations are 

([/l(xi,tl),/l(x2,t2)]+)/2 = - X2,tl - ^2^ + G^,(X2 - Xi,t2 - tif] 

+ ^GAj (0, of [GAj (Xl - X2, il - t2) + Ga, (x2 - XI, t2 - ii)] , 
o 

(3-28) 

([/2(xi,ti),/2(x2,i2)]+)/2 = -y[Gyi,(xi -X2,ti -t2)^ + GA,(x2 -xi,t2 -ii)^], (3-29) 

([/3(xi,ti),/3(x2,i2)]+)/2 = — [G,1,(X1 -X2,«l -t2) +G^,(X2 -Xi,t2 -tl)]. (3-30) 

We can see that these quantities are symmetric under the exchange of the arguments 
Xl -s-^' X2,ti ^ t2- These correlations play an important role in our investigation of 
the relation between the fluctuation force fields in POM and those in IFM. 
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The imaginary parts of the second-order correlations are 

([/l(xi,tl),/l(x2,t2)])/2 = ^([0[]>(xi,ti),,^g>(x2,t2)]), (3-31) 

([/2(xi,tl),/2(x2,t2)])/2 = ^{[Po>{^l,tl)Jly{^2,t2)]), (3-32) 

([/3(xi,ti),/3(x2,t2)])/2 = y([0O>(xi,tl),<^O>(x2,t2)]). (3-33) 

These quantities are related to the coefficients of the Langevin equation. The 
Langevin equation can be expressed in terms of the above quantities as follows: 

□0<(x,t) + m^^{^,t) + dVj)(x - y)^\{y,t) + " y){4>l>{y,t))4><iy,t) 

+i J'ds J d^yid^y2Sf^{^-y2) 

X < ^ ( [4>t> (yus), 4>l^ (y2 , t)] ) 0< (y2 , t) 
+{[fi{yi,s),fi{y2,t)])0<{y2,s-t,t)] 

7 ( [A (y 1 , s) , /2 (y2 , t)] ) (y 1 , s - i , t) , 4>< (y2 ,t)]+ 



1 



z ( [/s (y 1 > > /3 (y2 , t)] > (y 1 , s - 1 , t) , 4>l (y2 , t)] - 



+g([</'o>(yi,s),'Ao>(y2,t)]+)[(/>o<(yi,s -t),0o<(y2)][(^<(yi,s - t,i),(A<(y2,i)]+ 

+^(['^o>(yi''5).'^o>(y2'*) - {'}'l>{y2,t))]+)[(p<{yi, s - t,t),(j)^{y2,t)]ip<{y2, s - t,t) 

= J d'y{h{y,t) + f2{y,t)4>o<{y,t) + h{y,t)4^{y,t)}6f^{^-y). (3-34) 

In case of the nonlinear equation, it is in general not clear what kind of relation should 
be satisfied between the coefficients of the Langevin equation and the fluctuation 
force field, lij In the present case, the fluctuation force fields are related to the 
coefficients of the Langevin equation that are similar to those given by the 2nd f-d 
theorem for the linear Langevin equation. eJ^'tZP For this reason, we regard these 
relations as constituting the 2nd f-d theorem for the nonlinear equation, ill 



§4. Semiclassical Langevin equation in IFM 



In this section, we present the result of the application of IFM to (j)'^ theory. 
It is essentially due to the work of Greiner et al.u) (For a detailed derivation, see 

*-* Shibata and Hashitsume studied the nonlinear Langevin equation in the quantal spin system 
interacting with a heat reservoir and derive a nonlinear fluctuation-dissipation theorem. lZP 

**^ It is worth noting that some terms of order on the l.h.s. of Eq. (3-34) cannot be expressed 
in terms of fi. 
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Appendix p.) 

The c-number field (p is divided into two parts again by employing the momentum 
cutoff Aj as 



X, t) = -i^ J d3p{</>(p, t)9iAr - IpI) + 0(p, t)^(|p| - Aj)}e'P^ 



(2vr) 

i<(x,t) + (/>>(x,t) 



(4-1) 



Here, the field 0<(x, t) is that of soft modes and (py{x,t) that of hard modes. 
The action of this system is 



to 



><V> 



+ 



+ 5, 



int L</^<, </^>J, 



where 



S. 



<i»\-V<{» 

Sint[4> 



* 1 

4 



)]= d x-(a^(/><(»a^ </.<(» - m 0<(»), 



(4-2) 



(4-3) 



* d'^x^icPt + 403 + Q^l(j)l + 40<</.^ + (Pi). (4-4) 



to 



4! 



Now, as in case of POM, we assume that the initial density matrix of the total system 
is given by the direct product of the system and the environment density matrices: 



p{(l><,<P'<,4>>,(l>'>;to) = Psi(t><,(p'<;to) <^ Phi4>>,4>'>;'to)- 

Then, the influence functional is defined as 



(4-5) 



Dcp, 



<t>>f 



D(P' 



xexp{i(50> [(/)>] + Sint[(p<,(t)>\ -5'</.> W>\ - Sint[(t>'<A'>\)}Ph{(t)>i,4>'>6to), 

(4-6) 

where the fields (p^i, <p'^^ and (p^j correspond to the fields i?!)>(x, to), ';^>(x, to) and 
(/)>(x, t), respectively. Sip is called the "influence action". Because it is difficult to 
calculate the influence action exactly, we usually carry out a perturbative expansion. 
In this calculation, we keep terms through the second order. Then, we have 



Re5/i.[0<, (/>'<] = / d^x<PA(x)h^{G\\{^)\<p{x) 
to 



to 



+X' I d^xi(i*X20(ti - t2) \ ^0^(xi)ImG^^+(O)Im[G^V(2^i " X2f]Hxi] 
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+ l<pA{xiWixi) + ^0i(xi)]Rc[G^^+(xi - X2M{X2M\X2 + l<pAix2))]'j , 

(4-7) 

Im5/F[0<,^<] = d'^xid'^X2 |-^^^(xi)Im[G^!^(xi - X2f]4>A{x2) 



--(/)zi(xi)(/)(xi)Re[G/+(xi - X2)'^](l)Aix2)(t>{x2) 



+ \4>a{xiW{x^) + ^4>\{x^)]lHGi'+{xi - X2)]ct>A{x2)[4>\x2) + ^</'i(x2)]| , 

(4-8) 



where 



4>{x) = ^(</'<(x) + ct>'^{x)), <t>A{x) = 0<(x) - (t>'^{x). (4-9) 
Here, we have introduced the propagator 

-X2) = i(T[0>(xi)0>(a;2)]) 

= e{tl - t2)i{cl)>ixi)cP>{x2)) + e{t2 - tl)i(0>(x2)0>(xi)) 

= 0{ti - t2)GAj{^i - ^2,h - t2) + e{t2 - h)GAj{yi2 - y-i,t2 - h). 

(4-10) 

Furthermore, we introduce the following notation for simplicity: 
Re[G^^^+(xi - X2r] = [e{h - t2) + - t2)] 

X^[G^,(X1 - X2,tl - t2T + {-ITGaA^2 - Xl,t2 - tlf], 

(4.11) 

MGi'+{xr - X2T] = [e{h - 12) - {-iTe{h - 12)] 

Xj.[GaMi - X2,tl - t2T - {-ITGaA^2 - XI, t2 - tlT]. 

(4-12) 

Note that the influence action Sjp is complex in general. Our next task is to 
define a new action Sip which is purely real. It is called the "stochastic influence 
action" and is defined as 



j Z?6^6^&i'i[6]^'2[6]^3[6]e'^'^t'^<'*<'^'l- (4-13) 



Here, the stochastic weights Pi are given by 

Art 
to 



PilCi] = iV.exp d^xid''x2UxiK\xi - X2Mx2)^ , (4-14) 
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where 



1 



,2i-l 



2 



(4-15) 
(4-16) 
(4-17) 



Here, Ni is a normalization constant. The auxihary fields ^1, ^2 and ^3 have been in- 
troduced to eliminate the imaginary part of the influence action. Then, the stochastic 
influence action Sjp is given by 

(4-18) 

To prove Eq. (4-18), we have to assume the relation 

,-{i/2)x-Ax ^ [deti2nA)]-y' J Z^^exp [-\eA'H + ^X^c) ■ (4-19) 



(See Appendix C for details.) 

Finally, the stochastic effective action is defined as 



(4-20) 



Note that this action is real. If the soft modes behave quasiclassically, we can apply 
the variational relation 



5S, 



effl 



S4>A 



(4-21) 



to the stochastic effective action to obtain a semiclassical Langevin equation. In 
order to carry out this functional derivative, we must define a functional derivative 
for the soft modes. Greiner et al. used the usual definition, 



(x 



6^^\x'-x). 



S(j)<^{x 

Thus we obtain the semiclassical Langevin equation 

□0<(x,O + mV<(x,t)+ E (-l)^54%llM 

i=a,-,f ^ 



(4-22) 



[6 (X, t) + (X, t)^2 (X, t) + (l>l (X, t)^3 (X, t)] , 



(4-23) 



16 



T. Koide, M.Maruyama and F.Takagi 






(a) 



(b) 



(c) 



(d) 




(e) (f) 

Fig. 2 :Diagrams that contribute to the semiclassical Langevin equation in IFM. 
The dashed hues represent the soft modes and the sohd hues the hard modes. 



where 
(5Re[5 



Si 



6Re\S 



if\ 



64>a{x) 



5Re\S 



IFl 



64>a{x) 



6Re\S 



6(f)A{x) 



dReiS^^] 



^A{x 



A 



><(x,t), 



-yG^,(O,O)0<(x,i) 



(4-24) 
(4-25) 



t-to 



dr / dV— G^,(O,O)[G^,(x-x',r)2-G^^(x'-x,-r)2]0<(x,t), 
J 4 



(4-26) 



t~to 



dr / dV— [G^,(x-x',rf-GA,(x'-x,-r)3]0<(x',t-r)(4-27) 
J 

</.<(x) J^''\t J d3x'!^[G^^(x-x',r)2-G^^(x'-x,-T)2]</>2^(x',t-r), 

(4-28) 



t-to 



<Pi,{x) / dT / dV— -[G^,(x-x',r)-GA,(x'-x,-T)]0^<(x',t-T 



12 



(4-29) 



Here, the superscripts (a),(b),- • • correspond, respectively, to the diagrams (a),(b),- • • 
in Fig. 2. This classical treatment of the soft modes is justified in weakly coupled 
bosonic quantum fields at high temperature. (See Refs. ||), ^ and ^T|) for details.) 

In IFM, the derived semiclassical equation in general possesses a time-convolution 
integral. In the present case, it is contained in Eqs. (4-27), (4-28) and (4-29). To re- 
move such an integral, Greiner et al. introduced an approximation called the "linear 
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harmonic approximation" . To explain this approximation, we consider the equation 



di 



a{t) = —iu>a{t) + g dsF{s)a{t — s) 



(4-30) 



for some variable a{t), where w is a frequency, g is a coupling constant and F{s) is 
an appropriate function of s. This equation has a time-convolution integral in the 
second term on the r.h.s. The quasi-instantaneous approximation is often used to 
remove such an integral. In this approximation, one makes the replacement 

da{t) 



a{t 



a{t) 



dt 



(4-31) 



This approximation is justified when the time dependence of a(t) is sufficiently weak. 
0) However, if we use this approximation in t/)^ theory, the dissipation effect is lost. 
To retain the dissipation effect, an additional approximation is needed. For example, 
one may use a propagator with an explicit width instead of the propagator defined 
in Eq. (4-lO).0)'i) 

The linear harmonic approximation recently introduced by Greiner et al.& im- 
plies the replacement 

a{t - s) ^ ait)e''^' . (4-32) 

Here, the dependence of a{t — s) on the integration variable s is replaced with that 
of the free vibration. It is obvious that the time-convolution integral in Eq. (4-30) 
is eliminated through this replacement. In this approximation, we can retain the 
dissipation effect without imposing any further approximation in (p^ theory. However, 
the physical meaning of the approximation is not clear. 

In (p'^ theory, the linear harmonic approximation implies the replacement 



- (k, t) cos uj]^T — (/)< (k, t) — sin uji^t. (4-33) 



With this replacement, the three terms defined in Eqs. (4-27), (4-28) and (4-29) can 
be rewritten as follows: 



5Re[5^'^^ 



IF\ 



5(pAix) 



dT 



d3x'^[G^,(x',r)3 - G^,(-x', -r)3](^<(x - x', -r,t), (4-34) 



5(t)A{x) 



5Re[5}5] 



t-to 



dT 



d3x'^[G^,(x',r)2 - G^,(-x', -r)2]^<(x - x', -T,t)\ 



t-to 



(4-35) 

dr ( dV^[G^,(x',r)-G^,(-x',-T)](^<(x-x',-T,t)3. 



12 



(4-36) 



Here, we have introduced the function [cf. Eq. (3-18)] 



(^<(x,T,t) 



(27r)3 Jo 



1 



d k((/)<(k, t) cos WkT -|- 0<(k, tj — smtiJkTje 

Wk 



ikx 
(4-37) 
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Finally, the semiclassical Langevin equation without a time-convolution integral is 
obtained if we use the above expressions instead of Eqs. (4-27), (4-28) and (4-29). 

Now, we examine the r.h.s. of Eq. (4-23). In the IFM calculation, the auxiliary 
fields S,i are interpreted as the fluctuation force field.il''© Their expectation values 
are calculated by using the stochastic weights Pi. Then, the first order correlations 
are 



(ei(x,t))p = (6(x,t))p = {u^,t))p = 0, 



where 



(ei(x,t))j 



DCiPi[mii^,t). 



(4-38) 



(4-39) 



To know the characteristics of the correlations, it is necessary to calculate the cor- 
relations up to at least second-order. We have 



(^l(xi,ti)^l(x2,t2))p 
(6(xi,tl)6(x2,t2))p 

(6(xi,ii)?3(x2,i2))p = 

(Ci(xi,tl)^j(x2,t2))p = 0. 



— [Gyl^(xi - X2,il - ^2)^ + G'yl^(x2 - Xi,t2 " ^l) 

A2 



(4-40) 



[GAii^l - X2,tl - t2) + G^^(X2 - Xi,t2 - tl) ], 

(4-41) 

[G^,(xi -X2,ti -t2) +G'^,(x2 -xi,t2 -ti)](4-42) 



(fori/j) 



(4-43) 



It is noteworthy that there is no correlation between and for i ^ j- 

In §3, we have pointed out that, in POM, the 2nd f-d theorem is satisfied for 



at there is a 
9 To confirm 



the coefficients and the fluctuation force fields. In IFM, it is said 
relation called the "generalized fluctuation-dissipation theorem". c3) 
the validity of this relation, it is convenient to introduce the following functions 
7WW(k,cj): 



d^k 



M=0 

kx 



(27r)2 



e{Ai 



k|) I dTv.<(k,-T,t) j — A^W(k, 



UJ) 



(4-44) 



4><{x) 

ft-to 



J2^^ Jo (2^ 



k-kil) 



dr(^<(k-ki,-r,t)^<(ki,-T,t) j —e-^^^M^^\k,uj), 
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(4-45) 



<^<(^) / 7^:^^^''" ( " ^7nir (^i^i - Ik - ki - k2|) 



{2^^ Jo ~Wf 



(4-46) 

Here, y5<(k, r, t) is the Fourier transform of (p^{x,T,t). Furthermore, the Fourier 
transform of the function Ii{x), which is the inverse of the matrix I^^{x) given by 
Eq. (4-15), (4-16) and (4-17), is defined as 



Then, we obtain the relation 

Ii{k,uj) = -K»H7W«(k,u;), (4-48) 



id 

where 



K'^'' (uj) oc u;^—^ . (4-49) 

exp fiio — 1 

This relates the quantities which are the coefficients in the Langevin equation, 

with the quantities Jj, which characterizes the fluctuation force field. This is the 
generahzed fluctuation-dissipation theorem. In Refs. 42) and ^), it is pointed out 
that this relation is equivalent to the well-known Einstein formula in the classical or 
high-temperature limit for the case of quantum Brownian motion. 

We now discuss the arbitrariness of the definition of the fluctuation force field 
in IFM calculation. In the case of IFM, an auxiliary field is introduced to define 
the real stochastic effective action. It is interpreted as the fluctuation force field. 
However, in such a definition, this force field cannot be uniquely fixed. For example, 
the influence action does not change even if the term SI~^, which has the property 

6ir\x^-X2) = -6I-\x2-x^), (4-50) 



is added to /■ ^, because of the property 

rt 
'to 



f <fxi(fx2Uxi)5KHxi - X2)UX2) = 0. (4-51) 
J to 



Another point of arbitrariness in the definition of the fluctuation force field is 
that it is possible to use only one kind of auxiliary field instead of three. EP In this 
case, the stochastic influence action is given as 



SiF[4><A'<;i'\ = MSif[(P<A'<]] + t d''x{^^{x) - (t}'^{x))i'{x), (4-52) 



to 
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and the stochastic weight is 

P'i^] = Ne^py^'d^xid^X2^\xi)I-\xi - X2)^'{X2)^ , (4-53) 

where 

I~^{xi - X2) = Ii^ixi - X2) + ^{Xi)l2^{xi - X2)4>{X2) 

+ [4>\xi) + ^<l>\{xi)]I^\xi - X2)W{X2) + ^<I)\{X2)] 

= Ii^{xi - X2) + (f><:{xi)l2^{xi - X2)(t><{x2) 

+cl>l{xi)I^\xi - X2)4>l{x2). (4-54) 

In the last hne, we have taken the hmit (j)A{x) — 0. The derived semiclassical 
Langevin equation is 

□</><(x,t)+m20<(x,t) + ^4(x,t)+ E (,,^ ^Re[5g[^<,0-<]] 

i=a,b,c,d,e ^ <^^^P 

= e'(x,t). (4-55) 

Now, we have only one kind of fluctuation force field. Therefore, the fluctuation force 
field introduced in Eq. (4-23) is only one of several possibilities. There is still another 
point of arbitrariness, as is discussed in the next section. It plays an important role 
when we compare the quantum Langevin equation derived using POM with the 
semiclassical Langevin equation derived using IFM. 

§5. Quantum Langevin equation in POM and semiclassical Langevin 

equation in IFM 

To this point, we have applied POM and IFM to (f)^ theory and have derived 
two Langevin equations, (3-14) and (4-23). The equation derived with POM is 
the quantum Langevin equation, while that derived with IFM is the semiclassical 
Langevin equation. It is useful to compare the POM quantum Langevin equation 
with the IFM semiclassical Langevin equation and study the correspondence be- 
tween them. For the remainder of this section, we ignore the difference between 
operators and c-numbers for the soft modes. More precisely, we make the identifica- 
tion [A, = 2AB, where A and B are appropriate c-numbers corresponding to A 
and B, respectively. 

First, we compare the l.h.s. of the two equations. It is easily seen that there are 
two differences between the two equations. The first difference is that in POM, there 
are additional contributions that correspond to the diagrams shown in Figs.l (g) and 
(h), which contain the propagator of the soft modes. Such terms are not included in 
IFM equation, because the quantum effects from the soft modes are omitted there. 

The second difference is that between the cutoff delta function and the Dirac 
delta function. The origin of this difference is the definition of the functional deriva- 
tive, which is used to derive the Langevin equation from the stochastic effective 
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action in IFM. Greiner et al.cP used the usual definition of the functional derivative, 

^<P<{x') [</.<(x')+e<5(4)(x'-x)]-(^<(x') 
^, , s = lim 

= 5^^\x'-x). (5-1) 

In this definition, the function with small variation (/)< {x') + e^*^*^) (x' — x) includes 
the hard modes. This is unsatisfactory because, in our method of calculation, all the 
fields that contain the hard components should be integrated or projected out. This 
point can be confirmed by comparing the Fourier components of the two Langevin 
equations. We consider, for example, the third terms on the l.h.s. of the two equa- 
tions. In IFM, the k component of the third term is 

^ / d3xc/>3<(x)e-^l'- 

^ r'd3k^d3k2^(yl/-|k-ki-k2|)</.<(ki)0<(k2)</><(k-ki-k2), 
JO 



6(27r)f 



and in POM, it is 



(5-2) 



j d^x j dV<5S(x-y)</<3<(y)e-kx 

l"^' d^\iid^\i29{Ai - \k\)e{Ai - |k - ki 



6(27r)6 

X0<(ki)(/)<(k2)0<(k - ki - k2 



(5-3) 



Here, the time argument of (/>< has been omitted. The only diff'erence between 
Eq. (5-2) and Eq. (5-3) is the existence of the step function 9{Ai — |k|) in Eq. (5-3). 
The absolute values of the momenta ki , k2 and k — ki — k2 are always smaller than 
Aj in both cases. However, the situation is different in the case of the momentum 
k. In POM, only components smaller than Aj are included, even for momentum k. 
Contrastingly, in IFM, all components satisfying the condition < |k| < 3/1/ are 
included, because the step function 9{Ai — |k|) does not appear. This means that 
the separation of soft modes and hard modes is incomplete in the IFM calculation. 
To solve this problem, we use the following definition of the functional derivative: 

^^ = 5(x'o-xo)5S(x'-x), (5-4) 

where Jj^^ (x' — x) is the cutoff delta function defined in Eq. (3-8). If this new 
definition is adopted in IFM, the l.h.s. of the two equations coincide, except for the 
terms corresponding to (g) and (h) in Fig. 1. 

Now, we compare the r.h.s. of the two equations. If we use the definition (5-4) 
instead of Eq. (4-22), we obtain Eq. (4-23) with the r.h.s. replaced by 

/ dV{6(y,t) + <A<(y,i)6(y,i) + </'<(y,tK3(y,t)MS^(x-y). (5-5) 
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Comparing this with the r.h.s. of Eq. (3-14), we see that there are stiU the foUowing 
differences between them. First, the time dependences of the soft modes that appear 
on the r.h.s. of the two equations are different. In IFM, this time dependence is 
determined by solving the equation, while in POM, they evolve as free fields. This 
is a considerable difference. In POM, the r.h.s. of the quantum Langevin equation 
vanishes when we calculate the expectation value with the initial density matrix, 
because the soft modes cannot include the fi dependence. By contrast, in IFM, 
the soft modes included in the r.h.s. of the semiclassical Langevin equation have a 
dependence, and therefore the r.h.s. does not necessarily become zero when we 
calculate the expectation value with stochastic weights; that is, 

(/,(x,t)(^^o<(x,t)) =0, (5-6) 
(C,(x, t)r< (x, t))p = (^^(x, t)r< (x, t; Ci))p + 0. (5-7) 

Second, the correlations of the fluctuation force fields are different. Consider, for ex- 
ample, the second order correlations (^i(xi, ii)^j(x2, t2))p and (/i(xi, ii)/j(x2, i2))- 
The former is real and invariant under the exchange of the arguments xi -(-^ X2 , 
til as seen from Eqs. (4-40), (4-41) and (4-42). The latter is complex and is not 
invariant under the same exchange, as seen from Eqs. (3-22), (3-23) and (3-24). 
Therefore, the corresponding correlations in POM are considered to be the sym- 
metrized correlations given by the real parts of the second-order correlations defined 
in Eqs. (3-28), (3-29) and (3-30). Then, we notice that Eqs. (4-41) and (4-42) are 
identical to Eqs. (3-29) and (3-30), respectively: 

([/2(x,t),/2(x',t')]+)/2 = (6(x,t)6(x',t'))p, (5-8) 
([/3(x,t),/3(x',t')]+)/2 = (6(x,t)6(x',t'))p- (5-9) 

Contrastingly, Eq. (4-40) differs from Eq. (3-28) because of the existence of the 
additional term 81\ : 

([/i(x,t),/i(x',iO]+)/2 = (6(x,i)a(x',0)p-^'^^i(x-x',i- 0,(5-10) 

where 

(x - x', i - t') = -^X^A^ ^'^^^ (x - x', i - i') + Ga, (x - X, - t)\ 

= yG^^+(OflmGjV(x-x'). (5-11) 

However, this difference can be removed by modifying ^i. First, note that bl\ has 
the following property 

^ d'^Xid^X2(j)<ixi)SIi\xi - X2,tl - t2)(l><{x2) = 0. (5-12) 

This property can be proved using the fact that (t)<{x) contains only soft modes, 
while 51^^ contains only hard modes. With this fact, the imaginary part of the 
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influence action (4-8) can be rewritten as 
ImS'7F[0<,^<] 

^ J' d^xid^X2 l^ct^Aixi) (^-^Im[G^^+(xi - X2f] - lGi'+(,OflmG%{xi - ^2)) 
~(l)A{xi)^{xi)Re[G^'_^{xi - X2)'^](l)A{x2)4'ix2) 

+IMxi)[4>''{x,) + ^0i(xi)]Im[G^V(xi - X2)]Mx2)[4'''{x2) + ^<^i(x2)]} . (5-13) 

This imaginary part can be reproduced by replacing ^1 in the stochastic effective 
action (4-20) with the auxiliary field with the stochastic weight 

PiWi] = K exp d^xid^X2^[{xi)iir\xi - X2) + Sir\xi - X2))el(^2)| . 

(5-14) 

The second-order correlation of ^( then becomes 

= —ihixi - X2) + Sh{xi - X2)) 

= ([A(x,t),A(x',0]+)/2. (5-15) 

In short, we have shown that the second-order correlations in IFM coincide with 
those in POM by utilizing the arbitrariness of the definition of the fluctuation force 
field in IFM. However, it should be noted that in POM, there is also the nonzero 
correlation 

([/l(xi,il),/3(x2,i2)] + )/2 = -- ^G^,(0,0)[Ga,(xi - X2,il - t2) + G^, (x2 - Xi,t2 " tl)]- 

(5-16) 

The corresponding correlation is implicitly assumed to vanish in IFM. 

Now, we briefly discuss the fluctuation-dissipation theorem. In IFM, it is pointed 
out that the generalized f-d theorem is satisfied. This theorem asserts that there is 
a relation between the quantities A4^^\ which are defined by Eqs. (4-44), (4-45) 
and (4-46), and the quantities /j that determines the second-order correlation of the 
fluctuation force field. In order to see whether such a relation is satisfied in POM, 
it is necessary to identify the corresponding quantities M.^^^ and li in POM. From 
the fact that the l.h.s. of the two equations have the same form, it is seem that the 
in IFM are the same as those in POM. On the other hand, the correspondence 
between the r.h.s. of the two equations is not clear. We assume that /i, I2 and I3 
in POM correspond to the l.h.s. (or the r.h.s.) of Eqs. (3-28), (3-29) and (3-30), 
respectively. Equations (3-29) and (3-30) arc proportional to I2 and given in IFM. 
Therefore, the generalized f-d theorem is satisfied for them. However, the r.h.s. of 
Eq. (3-28) is not proportional to Ii, and therefore, the generalized f-d theorem is 
not satisfied for the fi part. 
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In POM, there is a relation (the 2nd f-d theorem) between the fluctuation force 
fields and the coefficients in the quantum Langevin equation. There is no such rela- 
tion in IFM, because the imaginary part of the second-order correlations is absent. 
In other words, the fluctuation force fleld in POM possesses more information than 
that in IFM. 

§6. Derivation of a semiclassical Langevin equation in POM 

In this section, we derive a semiclassical Langevin equation from the quantum 
Langevin equation (3-14). 

First, we must remove quantum effects that come from the system. For this 
purpose, we replace the operators of the soft modes with c-numbcrs. There arc two 
possible methods for this replacement. In the first, we make this replacement after 
taking the expectation value using the initial density matrix: 

^< (Xl , 1 1 ) <^< (X2 , 1 2 ) • • • (Xn , ) Tr [/9^< (xi , 1 1 ) (X2 , 1 2 ) " " " (Xn , ) ] 

~ <?i><(xi,tl)^<(x2,t2) •••'?!'< (Xn,tn), (6-1) 

for n = 1, 2, 3, • • •. Here p is the initial density matrix of the total system. In the 
second, we replace the operators themselves with the c-numbers: 

0< (Xl , ti)4>< {■X2,t2) ■ ■ ■ 4>< (Xn, tn) 0< (xi , ti)4>< {yi2M) ' ' ' <t>< (Xn, ^n) • 

(6-2) 

In the former method of replacement, we cannot obtain the fiuctuation force field, 
because the r.h.s. of the equation becomes zero when we take the expectation value, 
as follows from its definition. For this reason, we adopt the latter method of replace- 
ment. Note that a similar replacement should be made for 0o<(x, and ^o<(x, in 
the fluctuation force field. Furthermore, in Eq. (3-14), there are two terms including 
the soft mode propagator, which correspond to the diagrams shown in Figs. 1 (g) 
and (h). In deriving our semiclassical equation, we simply drop these terms. 

Next, we have to replace the fluctuation force field with a c-numbcr. The simplest 
method to do so is to introduce c- number fields /i , /2 and /s that have the following 
correlations: 

(/l(x,t)), = (/2(x,t))e = (/3(x,t))e = 0. (6-3) 
(/l(xi,tl)/l(x2,t2))c = -^{'^J^A'^^ - ^2, ti - tif + Gaj{^2 - Xl, t2 - hf] 

+^Gai (0, 0)2 [Ga, (xi - X2, ii - t2) + Gaj (x2 - Xl, i2 - ii)], (6-4) 
o 

\2 

(/2(xi,tl)/2(x2,t2))c = - — [GaA^I - X2,«l - t2f + Gyl, (X2 - Xi,t2 - tl)'],(6-5) 

(/3(xi,il)/3(x2,t2))c = -^[Gyl^(xi -X2,tl -12) + Gyl^ (x2 -Xi,t2 " ^l)], (6-6) 

\2 

(/l(xi, tl)/3(x2, t2))c = -^Gaj (0, Q)[Gaj (Xl - X2, h - t2) + Gaj {^2 - ^iM - h)]- 

(6-7) 
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Here, ( )c represents an average evaluated with respect to suitable stochastic weights. 
Other second-order correlations are assumed to be zero. These values are determined 
so as to reproduce the real part of the second-order correlations, in other word, the 
symmetrized correlations of /j. If there are no correlations between different /j, we 
obtain Gaussian stochastic weights, as in case of IFM. Actually, there is a correlation 
between f\ and /s in this case, and it is not clear whether or not one can construct 
a stochastic weight that reproduces all the above correlations. 

Because of the above stated problem, we consider another replacement, which 
is explicitly defined in terms of a definite stochastic weight. In the IFM calculation, 
the field is treated as the stochastic variable. We regard the creation-annihilation 
operator as the stochastic variable. The creation-annihilation operators of the hard 
modes, which develop as free fields, have the correlation 

(a|c(i)«k'(i')) = 'IV£[pBaj^ak,e^-k*e-^-k'*'] 
= n^^e^'^k(*-*')<5{3)(k_k/) 

= n(k,k';i,t') (6-8) 

We consider the replacement of the creation-annihilation operators of the hard modes 
ak(t) and a\j^t) with the c-numbers P\^{t) and respectively. Then, to realize the 

above correlation for these c-numbers, we introduce the stochastic weight P[/?, /?*] 
which has the property 

(/3kW/3k'(i'))^= n [D/3kD/3*^P[l3,l3*]l3*i,{t)P^,{t') 

= n(k,k';i,t'), (6-9) 



where 

P[P,P*] = A^exp 



- dsids2 J (fki(fk2l\^{si)n ^{ki,k2; si, S2)l3y^^{s2) 

(6-10) 



In this replacement, the field operator ^o>(x, t) is replaced with the c- number field 
(j)>{f3,(3*,x,t), given by 

</>>(/?, /3*,x,i) = /^d3k^=i==[/3k(t)e*- + /3i^(i)e-*-]. (6-11) 

Using this function we replace fi with c-numbers as follows: 

/i(x,t) /i(/3,/3*,x,t) = ^</.>(/3,/?*,x,tf , (6-12) 
/2(x,t) /2(/?,/3*,x,t) = ^(0>(/3,r,x,t)2 - (<^>(/3,/3*,x,i)2)^),(6-13) 
/3(x,t) /3(/3,/3*,x,t) = ^<^>(/3,/?*,x,t), (6-14) 
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Summarizing the above result, the semiclassical Langevin equation in POM is 
given by 



□0<(x,t)+m2</><(x,i)+ /tiVS(x-y)G4</><(y,0] 

i=a.,---J 

= J d'ysf^ (x - y){/i (/?, r , y, t) + <^o< (y, t)f2{P, (5*y, t) + 4>l< (y, t)h (/?, /5*y, t)}. 

(6-15) 

For comparison, the semiclassical equation in IFM is given by 

(x, t) + n^c^^ (x, t)+ J2 J '^'y^^'^ - [(/)< (y , i)] 

i=a,-,f 

= I dW)(x-y){Ci(y,t) + <^<(y,t)6(y,t) + 0|(y,t)6(y,i)}, 



(6-16) 



where 

Ga[My,t) 
Gb[<P<iy,t) 

Ge[My,t) 
Gf[My,t) 



-^G^,(0,0)</><(y,i), 

k2 rt 



(6-17) 

(6-18) 



--^ ds'l d^yi[GA^{yi-y,s-tf-GA,{y-yi,t-s)^]GAM(i)4><{y,t), 

(6-19) 

-^J^ds J d^yi[GA,{yi-y,-sf -GA,{y -yi,sf]ip<{yi,-s,t), (6-20) 



= — ^ y d^yi[GAi{yi - y,-sf - GAiiy - yi, sf]ipl.{yi,-s,t)(f)<{y,t), 



(6-21) 



A2 
12 



J^dsJ d^yi[GAj{yi - y,-s) - GAj{y - yi,s)](p%{yu-s,t)^%{y,t). 

(6-22) 



It is clear from the discussion in the previous section that the l.h.s. of the two 
equations arc identical except for the difference between the Dirac delta function 
and the cutoff delta function. For the fluctuation force fields, which are the r.h.s. 
of the equations, there is a difference between POM and IFM at the level of the 
semiclassical Langevin equation: In IFM, ^i, ^2 and ^3 are treated as stochastic 
variables, while in POM, there are only two stochastic variables /? and /?*. As a 
result, there are some differences in the correlations. To show this, we calculate the 
correlations of /j (/?,/?*, x, t). 

The first order correlations are 



(/i(/?,/?*,x,t))^ = (/2(/3,r,x,0)/3 = (/3(A/?*,x,0)/3 = 0. (6-23) 
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Thus all the first-order correlations disappear. This is also the case in IFM. The 
second order correlations are 

(/i(/3,/3*,xi,ii)/i(/3,r,X2,t2))/3 

1 



/ d\d%d% ' no;, no;, n^, 

3 Jai 8{2TTfoJi^uJi^uJi^ h h h 

X cosa;i^(ti - t2) cosu;i2(ti - ts) cosa;i3(ti - i2)e'(^i+^2+^3)(xi-x2) 

+2A2 d\dH2d\——^ n^, n^, n^, cosiUi^{ti - t2)e'^'^'''-''^\ 

JAj %{2'Kfuj\^LJl^LJl^ ll I2 is 3V 

(6-24) 

(/2(Ar,Xl,tl)/2(/3,/3*,X2,t2))/3 

= 2A^ / (j^lid^l2 ,,^ n^, n^, 

xcosiOi^iti-t2)cosoJi^iti-t2)e'^'^'+^^^^'''-''^\ (6-25) 

(/3(Ar,Xi,ti)/3(/3,/?*,X2,t2));3 

= n / ^m)/3(/?,r,Xl,tl)/3(/3,r,X2,i2) 

*'2(27r)3a;k"'"'^k 

(/l(/3,r,Xi,ti)/3(/3,r,X2,t2))/3 



= T / 2n^ cosa;k(ti - i2)e^''("^-"^), (6-26) 

(/3,r,Xi,ti)/3(/3,r,X2,t2))/3 

/ d3lid3l2— — n^, cosui^{ti-t2)cosiOi^{ti-t2)e'''^'''-'''\ 



(6-27) 

These correlations arc identical with the symmetrized second-order correlations if 
l + 2n^^ is replaced with 2n^.p in Eqs. (3-28), (3-29), (3-30) and (5-16). This can be 

understood as follows. The creation-annihilation operators and depends on 
their order of the operation: 

{aW') / (ak'ajc)- (6-28) 

However, in the replacement under consideration, the corresponding situation does 
not exist, as 

{PlMp = (/3k'/3k)/3- (6-29) 

In other words, in the present replacement, all the contributions from the vacuum 
are ignored. 

In this section, we have derived a scmiclassical Langevin equation from a quan- 
tum Langevin equation in POM ignoring the noncommutativity of the soft mode 
operators and diagrams with soft mode propagation. The problem we encounter is 
that the procedure of replacing the fluctuation force fields with c-number stochastic 
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variables is not unique. With the procedure we employed here, the correct correla- 
tions between the fluctuation force fields with different indices, like Eq. (6-27), are 
realized but the correct contribution from the vacuum is not realized. 

§7. Conclusions and discussion 

In this paper, we applied POM to 0^ theory and derived correct quantum and 
semiclassical Langevin equations that describe the time evolution of the soft modes. 
We have examined in detail the differences and similarities of the results obtained 
using POM and IFM. 

In the case of POM, the quantum Langevin equation is derived systematically 
from a pcrturbativc expansion of the Hciscnbcrg equation of motion. The fluctuation- 
dissipation theorem of second kind is applied to determine the order of the expansion. 
In this formalism, we can derive an equation of motion with no time-convolution 
integral term without imposing an approximation, like the linear harmonic approxi- 
mation. To derive the semiclassical Langevin equation from this equation of motion, 
the quantum effect arising from the soft modes must be removed by replacing the 
operators of the system, i.e. the soft modes, with c-numbers, and we ignored the 
contribution that includes the soft mode propagator. Furthermore, we replaced the 
fluctuation force field with a c-number by regarding the creation-annihilation oper- 
ators of the hard modes as stochastic variables. 

In IFM, the semiclassical Langevin eqiiation is derived by applying a variational 
principle to the stochastic effective action. This equation has time-convolution in- 
tegral terms, and the linear harmonic approximation is employed to remove them, 
while maintaining the dissipation effect. To derive the stochastic effective action, the 
imaginary part of the influence action was removed by introducing auxiliary fields. 
These fields arc interpreted as the fluctuation force fields. 

We compared the semiclassical Langevin equation obtained using IFM with the 
quantum Langevin equation obtained using POM. If we ignore the fiuctuation force 
field, the difference between the two equations comes only from the difference be- 
tween the Dirac delta function and the cutoff delta function. This difference results 
from the definition of the functional derivative that is used to derive the semiclas- 
sical Langevin equation from the stochastic effective action in IFM. In POM, the 
momentum component that is larger than the cutoff Aj is always projected out, due 
to this cutoff delta function. In IFM, by contrast, the separation of the momentum 
is incomplete. Note that this agreement of the respective of the l.h.s. is also due to 
the linear harmonic approximation employed in IFM. There is no such agreement 
if the quasi-instantaneous approximation is used. In this sense, it may be the case 
that POM provides the basis for the validity of the linear harmonic approximation. 

With regard to the fluctuation force field, the relation between the two ap- 
proaches is more complicated. The soft modes that are coupled to the fluctuation 
force field in POM develop as free fields, and therefore the fiuctuation force field 
term always disappears if the expectation value is calculated using the initial den- 
sity matrix. Contractingly, in IFM, the time dependence of the soft modes is to 
be determined by solving the equation. Therefore, the expectation values of the 
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fluctuation force terms in IFM do not vanish in general. 

Moreover, in IFM, there is an arbitrariness in the auxiUary fields. It is note- 
worthy that the correlations (^i(x, t)^i(x',t'))p and ([/i(x, t), /i(x', t')]_(_)/2 can be 
made identical by utilizing this arbitrariness. In short, we can conclude that in 
POM, we can determine the fluctuation force field uniquely, and the result suggests 
that the fluctuation force fleld in IFM may have to be modified. In IFM, it is neces- 
sary to introduce another principle to determine the fluctuation force fleld uniquely. 
Furthermore, it is implicitly assumed that there is no correlation between different 
fluctuation force fields in IFM. Our POM study indicates that such an assumption 
is not necessarily valid. 

In short, the derivation of the Langevin equation in POM seems to be more 
systematic and less ambiguous than that in IFM. Furthermore, the fact that the 
differences between the two approaches can be attributed to ambiguities in IFM 
seems to indicate the superiority of POM. However, such a claim cannot yet be 
made conclusively, and it is a future problem to study the deeper physical meaning 
of the differences between the two approaches. 

Greiner et al. have asserted that the system should thermalize with the approx- 
imations used in the present paper. However, there is no common understanding 
about what condition a nonlinear Langevin equation should satisfy to ensure the 
evolution of the system toward thermal equilibrium. The well-known fiuctuation- 
dissipation theorem may ensure thermal equilibrium for the linear case, but this is 
an open question for the nonlinear case. c3 Therefore, it is not clear at present under 
what condition the equation of motion obtained with either POM or IFM guaran- 
tees the realization of thermal equilibrium. In the present paper, we have compared 
the two equations obtained from POM and IFM with respect to the assumptions 
and arbitrariness contained in each formalism. On the basis of this comparison, we 
conclude that POM is better than IFM. 

Finally, we derived the semiclassical Langevin equation in POM. A crucial prob- 
lem here is to determine how to replace the fiuctuation force field, which consists of 
an environment operator, with a c-number. In this paper, we replace the creation- 
annihilation operators of the environment with corresponding to the stochastic vari- 
ables. In this replacement, all the vacuum contributions are dropped. It is a future 
problem to further study the justification of this replacement. 

Appendix A 

Fluctuation- dissipation theorem of second kind 

In conventional statistical mechanics, the following two relations are known as 
the fiuctuation-dissipation theorem. The first is that the transport coefficient is 
represented in terms of the fluctuations of macroscopic variables. This is called 
the fluctuation-dissipation theorem of flrst kind (1st f-d theorem). The second is the 
fluctuation-dissipation theorem of second kind (2nd f-d theorem), according to which 
the transport coefficients are represented by the fiuctuation force. In POM, the 2nd 
f-d theorem is applied to determine the order of the expansion of the fluctuation 
force. 
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le 2nd f-d theorem is usually understood in terms of the Mori-Langevin equa- 



di' 



A{t) = iujA{t) - r dsifit - s)A{s) + f{t), (A-1) 
Jo 



where j4 is a macroscopic variable, <p{t) the memory function and f{t) the fluctuation 
force. The memory function is given by 

^(t) = (/(t),/t(0)), (A.2) 

where ( ) represents an inner product satisfying the following properties: 

(F,G)* = (G,F), (A-3) 
(^qF„G) = ^c,(F„G) (A-4) 

i i 

It can be seen in Eq. (A-2) that the memory function is represented by the time 
correlation of the fluctuation forces. This is the 2nd f-d theorem. From the above 
discussion, we can see that if the memory function is of second order in the coupling 
constant, the fluctuation force must be of first order for the 2nd f-d theorem to 
be satisfied. The Mori equation is derived exactly from the Heisenberg equation 
of motion, and therefore, the 2nd f-d theorem holds exactly in the case of a linear 
equation. In this paper, we would like to propose that this theorem holds to every 
order in the perturbative expansion in the case of a linear equation. 

In the calculation of the Langevin equation in POM, we expand the factor 
gjLQ(t-to)__i___^ zeroth order in Lj and determine the form of the fluctu- 

ation force. 



QV{t,to)Q 

1 



giLo(t-to) 



1 + (C(t, to) - l)Q + mt, to) - 1) + {C{t, to) - l)QiV{t, to) - 
~ Qe*^o(*-*o)iLO(to). (A-5) 

Now, we confirm that this fluctuation force satisfies the 2nd f-d theorem em- 
ploying the following model. The Hamiltonian is 

+ / diOiob\uj)b{Lo) + g / dio{a^b{uj) + ab^ (lo)), (A-6) 



where both a and 6(a;) are boson operators, which satisfy the following commutation 
relations: 

[a,at] = l, [b{uj),bHu')]=6{LO-Lo'). (A-7) 

We consider the case in which a is the system operator and b{oj) the environment 
operator and divide the Hamiltonian into three parts: 

= LVaa^a, (A-8) 
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He = J du>ujb\Lo)b{uj) , 
Hi = g J dw{a^b{u;) + ab^u)). 
The projection operator is defined as 

PO = b{0\O\0)b, 

where 

6H|0)6 = 0. 

Substituting these ingredients into Eq. (2-19), we have 

^a\t) = iuaa^t) - 5' / rfa; "'^"""''^'~/ at(t) + f{t), (A-13) 

where 

f{t) =ig j dioe'^^b\ij). (A- 14) 

This equation can be rewritten as 

^at(i) = iu;„at(t) + ^*ds(0|[/(i),/t(,)]|o)e--'^(*-^)at(t) + /(i) 
ft 

= iu;aa\t)+ / ds(t){t-s)e-''^''^^-'^a\t) + f{t) 
Jo 

= iujaa\t)+ I ds4){s)e-^'^'a\t) + f{t). (A-15) 
Jo 

This equation is the TCL equation, and therefore, we cannot compare this result with 

the Mori-Langevin equation, which has a time-convohition integral term. However, if 
we accept the approximation a\t — s) ^ e~'^^'^^a\t), as in case of the linear harmonic 
approximation, we can conclude that this equation satisfies the 2nd f-d theorem to 
this order of the perturbative expansion. Therefore, the definition of the fluctuation 
force (A-5) is consistent with the 2nd f-d theorem in this case of a linear Langevin 
equation. We apply this approximate fluctuation force to the nonlinear case. 



(A-9) 
(A-10) 

(A-ll) 
(A-12) 



Appendix B 

Derivation of a semiclassical equation in IFM 

For simplicity, we consider the case of quantum mechanics, in which the system 
is denoted by X and the environment by Q. The time evolution of the density matrix 
is described by 

PxvjQ{x,q\x\q';t) 

= J dxidqidx'idqp {x, q; Xi, qf, t, ti)pxuQ{xi, qf, x\, g-; ti)U*{x', q'; x-, g-; t, ti), 

(B-1) 
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where 



with 



U{x, q; Xi, Qi- 1, U) = r Dx' [' Dg'e^^(^''«')/?i, (6-2) 

Jxi Jqi 



S{x,q)=[ dsC{x{s),q{s)). 
Jti 

Here, the action S{x,q) is given by 

S{x, q) = Sx{x) + SQ{q) + Sint{x, q), 



(B-3) 



(B-4) 



where Sx{x) is the action of the system, SQ{q) the action of the environment, and 
Sint{x, q) the action of the interaction part. As in case of POM, we assume no initial 
correlation between X and Q: 

PxuQ{xi,qi;x'i,ql;ti) = px{xi,Xi;ti) ^ pQ{qi,q-;ti). (B-5) 

The reduced density matrix is defined by tracing out the environment degrees 
of freedom, as 

Pr{x,x';t) = Tr(q){pxvjQ{x,q;x' ,q';t)] 

J Jxi Jx[ 

dxidx'i Dx Dx'e'^-ff^'''''''>px{xi,x'i;t), (B-6) 

Jxi Jx'^ 

where the functions e^^"'^^'^"> and Seff{x,x') are the influence functional and the 
effective action, respectively. It is difficult to trace out the environment degrees of 
freedom completely, and therefore, we have carried out a pcrturbativc expansion 
in going from the first line to the second line. As an example, we consider the 
interaction 

Cintix,q) = xS{q). (B-7) 
The influence action Sif{x,x') is expanded up to second order in powers of x: 
Sif[x,x'] 

= dsF{s)[x{s)-x'{s)] 
Jti 



1 ft 

+- / dsds' 

2 Jti 



R{s, s') + il{s, s') ] x{s)x{s') + ( - R{s, s') + il{s, s) x'{s)x'{s'^ 



+ (^R{s, s') - il{s, s')j x{s)x'{s') +(^-R{s, s') - il{s, s')j x'{s)x{s') 



(B- 
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where F{s), R{s,s') and I{s,s') are real (Hermitian) functions. The effective action 
Sgff is a complex function. To make it real, we introduce a new variable ^, and with 
it we obtain the stochastic effective action, 

Seff{x,x',0 = Sx{x) - Sx{x')+Re{SiF{x,x')) + f' ds^{s)[x{s) - x' (s)] 

Jti 

= Sx{x)-Sx(x') + Sif(x,x',0- (B-9) 
This is derived by using the relation 

^i/hSeff{x,x') ^ j ]J^p^^-jQi/nS^ff{x,x',i) ^ (B-10) 

where P[$] is the stochastic weight defined by 

m = ^exp (^-^£dsids2asi)I-\s,,S2)as2)^ . (B-ll) 

In deriving this relation, we have assumed that it is possible to carry out the Gaussian 
integral. (See Appendix C for more details.) 

If the system behaves quasi-classically, the reduced density matrix becomes 
nearly diagonal. Then, the major contributing path is obtained by applying a vari- 
ational principle to the stochastic effective action: 



5Seff{x,A,0 



6A 

where 



= 0, (B-12) 

A=0 



X = , (B-13) 

A = x-x'. (B-14) 

The equation derived in this manner is a semiclassical equation that includes quan- 
tum effects of the environment: 



SSxjx) 
dx{s) 



+ F{s) + £ ds'R{s, s')x{s') = -C{s). (B-15) 



Here, we assume that the new variable ^ can be interpreted as a fluctuation force 
that is distributed with the stochastic weight P[^] defined in Eq. (B-ll). Then, 
Eq. (B-15) can be regarded as a kind of Langevin equation. 

Appendix C 

Gaussian integral 



To confirm Eq. (B-10), we calculate the following integral: 
J D^P[^]e^'^/^^^'^^^'^''^^ 
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Here, we omit the summation symbol. We assume that the matrix can be 
diagonalized by a unitary matrix U : 



(C-2) 



/ Ai \ 

V ■■• / 

The terms that appear in the brackets on the r.h.s. of Eq. (C-1) can be transformed 
as 

^C^ix^ - X[) - = -^vf + - yi) 



~{{m-^^{m-y'i)f + ^im-y'i}'}, 

(C-3) 



where {Ux).i = yi, (U^)i = rji. Furthermore, wc assume that all the eigenvalues Aj 
are positive. Then, we can carry out the Gaussian integral: 

I I' i__ _ /\\2-\ TT / I ^« 2 

'Vi 



2n 



[det{2hTrIij)]^/^ . (C-4) 



It follows that 



I £>^p[^]e(V't)5.Ha^y,^) = exp |^Re57F(x,x') - ^ {y^ - y'^ 

= exp ^ i^eSiF{x,x') + ^(xj - x'i)Iij{xj - 



where 



= [Det{2n7rl)]^/^. (C-6) 



Thus, the relation (B-10) is confirmed. Equation (4-19) is the quantum field theo- 
retical version of this relation. 



Projection Operator Approach to Langevin Equations in cp'^ Theory 
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